Evidence for Two— Band Superconductivity from Break Junction Tunneling on MgB 2 
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Superconductor-insulator-superconductor tunnel junctions have been fabricated on MgB2 that 
display Josephson and quasiparticle currents. These junctions exhibit a gap magnitude, A ~ 2.5 
meV, that is considerably smaller than the BCS value, but which clearly and reproducibly closes 
near the bulk T c . In conjunction with fits of the conductance spectra, these results are interpreted 
as direct evidence of two-band superconductivity. 

PACS numbers: 73.40. Gk, 74.50. +r, 74.70.Ad, 74.80.Fp 



The discovery of superconductivity in MgB2 has led 
to intense research activity, but the nature of the en- 
ergy gap, A, has been elusive. Tunneling spectroscopy, 
which is the most direct measure of this quantity, has 
revealed a large spread of A values and considerable 
variation in its spectral shape. Sharp, BCS-like tun- 
neling spectra have been observed in scanning tunnel- 
ing microscopy (STM) with a surprisingly small A = 2.0 
meV |ij . Other STM and point-contact studies revealed 
double-peaked spectra at low temperatures ]|,[| that 
were interpreted as evidence for two-gap superconduc- 
tivity. A provocative suggestion is that multiple gaps are 
a consequence of the coupling of distinct electronic bands 
j|. Our ability to fabricate superconductor-insulator- 
superconductor (SIS) break junctions has led to unique 
observations and we have gone beyond these initial re- 
ports to present more compelling evidence that MgB2 is 
one of the rare examples of two-band superconductiv- 
ity. In addition, our identification of a weak higher-bias 
spectral feature has provided important insight into the 
nature of the inter-band coupling. 

The simplicity of the crystalline structure in MgB2 al- 
lows for ab-initio calculations of its electronic properties 
, from which it is known that the Fermi surface consists 
of four sheets, two being two dimensional (2d) bonding a - 
bands and two being three dimensional (3d) bonding and 
antibonding 7r-bands. An and Pickett propose super- 
conductivity to be driven by the 2d cx-bands, where elec- 
trons are strongly coupled primarily to the E^g phonon 
mode. This raises important questions of how supercon- 
ductivity would manifest itself on the 3d sheets and how 
the tunneling density of states (DOS) would depend on 
the crystallographic orientation. 

A more recent work [|J treated the problem by reduc- 
ing it to two distinct bands which, in the clean limit, 
leads to the appearance of two isotropic gaps, A2 ~ 7.2 
meV and Ai ~ 2.4 meV, associated with the 2d and 
3d bands respectively. The small gap, Ai, on the 3d 
sheets is enhanced above its intrinsic value due to virtual 
phonon exchange (pair transfer) with the 2d sheets and 
should persist up to the bulk T c . The results of our tun- 
neling study address these issues in the following ways. 
First, the small gap feature is unambiguously tracked to 



high temperatures where it is still visible in the raw data, 
a key observation supporting two-band superconductiv- 
ity. These junctions only probe the band with the small, 
induced gap suggesting the SIS configuration strongly fa- 
vors tunneling between the 3d sheets. Second, a subtle 
spectral feature is observed in the conductance near 9 
meV that resembles strong-coupling effects. Using a the- 
oretical, two-band model |5j that treats pair transfer and 
the quasiparticle self-energy on an equal footing, we have 
quantitatively fit this feature. This indicates that self- 
energy effects originating from quasiparticle scattering 
between bands are important. Finally and importantly, 
using unique features of SIS junctions, we have more rig- 
orously ruled out proximity effects which otherwise can 
mimic the temperature dependencies and the quasiparti- 
cle self-energy effects of two-band superconductivity. 

Compact samples of MgB2 were formed from amor- 
phous B powder (4N's purity) and high purity Mg. The 
B powder was pressed into pellets under 6 kbar pressure. 
These free standing pellets were reacted with Mg vapor 
at 850°C for 2 hr in a BN container under 50 bar of 
Ar. During the diffusion reaction the pellets broke up 
into irregularly shaped pieces several mm on a side. The 
material typically showed T c — 39 K. To obtain a clean 
and smooth surface, the samples were polished until a 
shiny surface was exposed. No solvent was used and the 
samples were only cleaned in a flow of dry N2 gas. 

The tunneling measurements were performed on two 
different samples (A and B) using a point-contact ap- 
paratus |^] with a gold tip. This technique yielded 
superconductor-insulator-normal metal (SIN) junctions 
as well as SIS break junctions, the latter being confirmed 
by fits to the temperature dependent conductance as 
well as the observation of Josephson currents. Fig. [j] 
gives an example of a Josephson tunnel junction obtained 
on sample A. In the current voltage characteristic there 
is a distinct jump visible between a well-pronounced, 
slightly hysteretic Josephson current and the quasipar- 
ticle branch which displays a gap feature at eV — 2 A. 
In the M? vs. V spectra, sharp peaks are seen which can 
be fit to a BCS s-wave model with A = 2.1 meV and 
r = 0.5 meV. The fit is very good except for a weak 
spectral feature in the data near 9 mV, which we discuss 
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FIG. 1. Josephson tunnel junction on MgB2 displaying an 
energy gap of 2.1 meV. Upper inset: s-wave fit to the con- 
ductance characteristic. Lower inset: Comparison to an SIN 
junction on sample B showing a comparable gap magnitude 
at 1.6 K. 

later. The simultaneous observation of a well developed 
gap feature and a significant Josephson current estab- 
lishes these junctions to be of SIS geometry. The gap 
value is significantly less than the BCS expectation of 5.9 
meV, but consistent with the I c Rn product of 2.2 meV 
for this junction (other junctions gave comparable val- 
ues). Note that similarly low values of I c Rn were found 
in a recent study of Josephson junctions in MgB2 || pre- 
pared by a completely different break junction method. 
Although the focus of this study is on the SIS junctions, 
similar gap values were consistently found in SIN junc- 
tions (lower inset of Fig. |l|). Fits to a smeared BCS model 
yield A = 2.6 meV and T = 0.5 meV, that are consistent 
with the parameters obtained for SIS junctions. 

Many SIS junctions were obtained at different loca- 
tions on samples A and B. With increasing resistance 
of the junction the Josephson currents decrease, and for 
contact resistances higher than ~ 10 kfi no supercurrent 
is visible at zero bias. Nevertheless, the junctions could 
be easily identified as SIS by the evolution of the conduc- 
tance spectra with increasing temperature. SIN junctions 
show a rapid decrease in peak height at eV — A with in- 
creasing temperature, while the peak position remains 
approximately constant. In contrast, the gap structure 
of SIS junctions remains sharp even for elevated temper- 
atures, and the 2A-peak position follows the closing of 
the gap at T c . We traced three such junctions formed 
on different parts of the MgB2 sample A up to tempera- 
tures of ~ 30 K, and two of these sets of data are shown 
in Fig. ||, while the third is shown as a color map in Fig. 
[|. What is immediately evident in the raw data is that 
the small gap persists up to 30 K. Each curve shown in 
Fig. ^ was divided by Rn (the high bias resistance of the 
4.2 K data) and then fit to an s-wave gap SIS model us- 
ing the measured temperature. The smearing parameter 
r was adjusted once to reproduce the low-temperature 
data and then held fixed. This was done to prove the 
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FIG. 2. Temperature dependence of two high-resistance 
SIS tunnel junctions together with s-wave fits reproducing all 
major features of the conductance characteristics, including 
the evolution of a zero-bias peak due to thermally activated 
current. Each spectrum is rescaled by the indicated Rn- 

ability of the model to account for all major features of 
the temperature dependent data using A as the only ad- 
justable parameter. The fits shown in Fig. |^ reproduce 
the peak heights, the shape of the in-gap conductance as 
well as the evolution of the zero-bias peak that is due to 
thermally activated quasiparticle current. This zero-bias 
feature cannot appear in SIN tunneling but is a known 
feature of SIS tunneling. The magnitude of this central 
peak is surprisingly high, and this is caused by the un- 
usual case of a small gap which persists to comparatively 
high temperatures (much higher than the BCS T c con- 
nected with this gap) thus giving rise to unusually large 
thermal activation of quasiparticles. The BCS T c for a 
low temperature gap of ~ 2.5 meV is below 17 K, and the 
gap seen here is still clearly developed at a temperature 
around 30 K. This rules out a lower T c on the sample 
surface as an explanation of the small gap value as this 
would display a second transition near 17 K, which is not 
observed (see Fig. ||). Instead, these features indicate 
a bulk property is being measured that clearly deviates 
from strong coupling or BCS weak coupling theory. 

After taking the first set of temperature-dependent 
data (junction #1) the point contact tip was twice me- 
chanically retracted over ~ 100 /im and two new junc- 
tions (#2 and #3) were formed on the same sample (A), 
showing the identical spectral shapes and temperature 
dependence. Since the thread mechanism of the tip ap- 
proach does not preserve the microscopic lateral position 
on the sample, these junctions have to be regarded as en- 
tirely independent and their consistency thus proves the 
reproducibility of these observations. In addition it will 
be shown that SIS junctions on sample B display nearly 
identical low-T characteristics. Further tests, which in- 
cluded etching of the Au tip, suggested these junctions 
were formed by breaking off an MgB2 crystal fragment 
which then forms an SIS junction with the bulk material. 
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FIG. 3. Top panel: Color map of SIS junction conduc- 
tance spectra as a function of temperature. Bottom panel: 
Temperature dependence of the small superconducting gap 
Ai. Symbols are data extracted from SIS fits to the con- 
ductance spectra (different symbols represent three different 
experiments). The solid line gives A(T) according to BCS 
theory, using the expected ratio of 2A(0)/fcsT c . The dashed 
line gives a rescaled BCS dependence with both A and T c 
adjusted to fit the experimental data. 

This procedure of creating SIS break junctions is well 
known for point-contact tunneling into BSCCO Jl(]] . The 
ease with which reproducible SIS junctions were formed 
was quite surprising and this requires us to critically re- 
view our previously published data which showed similar 
spectral shapes at 4.2 K fli]] , but no T-dependence was 
measured. The data presented there were interpreted in 
terms of SIN tunneling with a gap value of 4.3 — 4.6 meV 
and no smearing parameter T. This analysis is similar 
in shape to an SIS characteristic with A = 2.5 meV and 
T = 0.5 meV (typical parameters used here to fit the 
data e.g. in Fig. g), leading to the possibility that some 
of the junctions in that work were in fact SIS type. This 
would bring the measured gaps of that work in line with 
what was consistently observed here. 

The values for Ai (T) extracted from the fits to the con- 
ductance data on these three junctions are shown in Fig. 
^. The uncertainty in the measured gap values is given 
by the error bars and there is excellent reproducibility 
among the three junctions. The solid line gives the BCS 
dependence for a low temperature gap of 2.5 meV with 
the expected T c = Ai(0)/(1.76fc B ) = 16.5 K. Compari- 
son of the data with the BCS fit shows a clear anomaly, 
the small gap persists up to temperatures far beyond 
the expected T c . The dashed line gives a rescaled BCS 
A(T) with the T c adjusted to fit the experimental data. 
This leads to an extrapolated T c near the bulk value for 
MgB2 but we want to strongly emphasize, that there is 
no basis for this type of rescaling within BCS theory. 
This should be contrasted with strong-coupling effects 
which commonly lead to a zero temperature gap that ex- 
ceeds the BCS value resulting in an enhanced gap ratio 
2A(0)/fcsT c . Here, this ratio is much smaller than the 



expected BCS value. 

We propose that the data of Fig. || together with our 
analysis of the conductance spectra presented below are 
compelling evidence of two-band superconductivity as 
suggested by Liu et al. Qj. The absence of any evidence 
for a second transition in Fig. [| indicates that Ai is pri- 
marily induced via coupling to the 2d band. We believe 
that the persistent and clear observation of only Ai in 
our SIS and SIN junctions is due to the much higher 
probability of tunneling into the 3d band. This band ac- 
counts for 58% of the total DOS Q which is insufficient 
to account for its preference. It rather is the dominance 
of tunneling by electrons with momenta normal to the 
barrier which favors the 3d band. Assuming a random 
orientation, there is a relatively low probability of being 
properly aligned with the 2d band. 

The low temperature conductance spectra reveal a 
weak but reproducible structure near 9 mV which we 
believe is related to the large gap, A 2 , but is not due 
to direct tunneling into the 2d band. This conclusion 
is supported by a calculation of the quasiparticle DOS 
on each sheet using the McMillan tunneling model ||. 
This model simulates a coupled, two-band system by 
including the BCS-type, virtual phonon coupling (pair 
exchange) between bands Jl2|,|l3| and also self-energy ef- 
fects from interband quasiparticle exchange. The model 
requires the solution of two simultaneous equations: 
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for the energy dependence of the two gaps, where the sec- 
ond equation is obtained by interchanging the subscripts 
1 and 2. These functions Ai^{E) are subsequently used 
to compute the DOS in both bands via the usual BCS 
expression. Convolution of these DOS then yields the 
desired SIS conductance characteristics. The model in- 
cludes six parameters: the intrinsic pairing gaps on both 
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bands in the absence of any interband coupling, A 
two scattering rates, I\2, related inversely to the times 
spent in each band prior to scattering to the other, and 
two smearing parameters, 2 , which were added to ac- 
count for lifetime effects within each band. Rather than 
treat each parameter as free, we fix the intrinsic pairing 
gaps by considering first principle calculations for MgB 2 
A? h 



viz. 



and AP h = 7.2 meV fL4). The first as- 



sumption is further justified by noting that only a single 
transition is observed in Fig. ^. The observed induced 
gap magnitude is then adjusted by appropriate choice of 
Ti : 2- These parameters are highly interdependent and 
cannot separately be determined. Good fits can be ob- 
tained for ratios T2/T1 between zero and ~ 0.5. We 
choose 0.25 as representative [ p~5| . The smearing param- 
eters finally arc needed to account for the broadening and 
the in-gap current. 

The low-temperature spectra from different junctions 
on samples A and B are shown in Fig. ^ to demonstrate 
the reproducibility of the data. Note that all 
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FIG. 4. Low temperature SIS tunneling spectra. The mea- 
sured data (dots) is normalized to a smooth background. For 
one set of data the fit (thin solid line) to a small (intrin- 
sic) BCS s-wave gap is shown along with a fit (thick solid 
line) to an induced gap. The parameters used are Aj h = 0, 
Af 1 = 7.2 meV, Ti = 4.0 meV, F 2 = 1.0 meV, V* = 0.5 meV 
and T2 = 1.0 meV. See text for details. 



rule out proximity effects. 

In conclusion, we have reproducibly observed a small 
energy gap, Ai ~ 2.5 meV, which smoothly closes near 
the bulk T c . This observation, along with a detailed anal- 
ysis of the conductance spectra, is indicative that MgB2 
is an example of the rarely observed phenomenon of two- 
band superconductivity. While the measured gap values 
are consistent with the first principles two-band model 
that assumes pair exchange between bands Q , our spec- 
tra are providing strong evidence that interband quasi- 
particle exchange is important. The data are inconsis- 
tent with surface layers of reduced T c or with proximity 
effects. 
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This research is supported by the U.S. Department of 
Energy, Basic Energy Sciences — Materials Sciences, un- 
der contract # W-31-109-ENG-38. 



characteristic features of the data taken on sample B, 
such as the gap magnitude and the additional structure 
outside the conductance peaks, coincide with those from 
sample A. For one junction (# 3), the result of a two- 
band fit is given along with the BCS model including 
a smearing of T = 0.5 meV. The BCS fit convincingly 
reproduces the in-gap conductance and the quasiparti- 
cle peak, but it cannot account for the features near 9 
mV, which are consistently seen in these SIS junctions. 
The two-band fit is close to the BCS behavior for low 
bias, but it also reproduces the shoulder and dip at the 
observed energy in agreement with the data. The higher 
energy spectral feature is an effect of the large gap, A2, on 
the quasiparticle self-energy in the 3d band DOS. Note 
that the dip drops below unity, a characteristic feature of 
quasiparticle self-energy effects, similar to phonon struc- 
tures, but which cannot be achieved by arbitrarily adding 
two BCS DOS from the two bands. The subtlety of the 
feature near 9 meV, as well as the dip, is more easily un- 
derstood as an intrinsic feature of the 3d DOS. Our data 
show no direct contribution from the 2d band |l6| . 

The same formalism used above to calculate the con- 
ductances is valid for a proximity sandwich, P, consisting 
of N and S layers that are coupled by a tunnel barrier 1 1| , 
i.e. P = NIS. Well-coupled N and S layers result in spec- 
tral features inconsistent with our observations [^7) . For 
P = NIS, the small gap in the N region has been shown 
to decrease exponentially with barrier thickness [Q. To 
explain the small gap of MgB2 as a surface proximity ef- 
fect (i.e., PIP' for our SIS geometry) would produce a 
noticeable splitting of the zero-bias peak (Fig. |J) unless 
the gaps in P and P' differed by less than ~ 20%. This is 
highly unlikely as it would require [ p"8| the tunneling bar- 
rier thicknesses in P for all MgB2 samples to be within 
~ 0.5 A. Thus it is our SIS geometry that allows us to 
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